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Abstract 

The 5 = 1/2 Heisenberg antiferromagnet in the ladder geometry is studied 
as a model for the spin degrees of freedom of SrCu203. The susceptibihty and 
the spin echo decay rate are calculated using a quantum Monte Carlo tech- 
nique, and the spin-lattice relaxation rate is obtained by maximum entropy 
analytic continuation of imaginary time correlation functions. All calculated 
quantities are in reasonable agreement with experimental results for SrCu203 
if the exchange coupling J ~ 850K, i.e. significantly smaller than in high-Tc 
cuprates. 
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The Cu-0 layers of SrCu203 have an internal structure of parallel double chains 
(ladders) .SH Cu spins within a ladder are exchange coupled with a strength expected to 
be comparable to that of high-Tc cuprates, whereas the inter-ladder coupling is weak, aris- 
ing from 90° Cu-O-Cu bonds. The spin degrees of freedom should therefore be well described 
by the Heisenberg model on a single ladder ,i defined by the hamiltonian 

i a=l,2 i 

where Sa,i is a spin-1/2 operator at site i of chain a. It is now well established that this 
system has a gap between the ground state and the lowest excitation for any ratio J2/ Ji 7^ 0. 
For J^ = J^ = J, the gap is A = 0.504J.§ 

Recent experiments on SrCu203 have been carried out by Azuma et a/.i and Ishida et a/..i 
Their results for the spin susceptibility x ^"^^ the ^^Cu NMR spin-lattice relaxation rate 1/Ti 
show clear evidence of a gap. Accordingly, the spin-echo decay I/T2G rate saturates at low 
temperatures, indicating a finite correlation length in the ground state.i However, comparing 
the data for x ^"^^ V^i with theoretical low-temperature results for the Heisenberg ladder 
obtained by Troyer et a/.,i there is a significant discrepancy; x indicates a gap A 420i^, 
whereas the behavior of 1/Ti suggests a gap close to TOOK.iili At first sight, one would tend 
to believe that the gap extracted from 1/Ti is the correct one, since the corresponding value 
of J ^ 2A is then close to the exchange constants typically found in planar cuprates. 

We have carried out quantum Monte Carlo (QMC) simulations of the Heisenberg ladder, 
and obtained results for the quantities discussed above. Here we present comparisons with 
the experimental results, and discuss a possible reason for the gap-size discrepancy found in 
earlier work. We argue that the formula used to extract the gap from 1/Ti is not applicable 
in the temperature regime where it was used, and that the gap obtained from x is more 
accurate. The calculated I/T2G is also in close agreement with the experimental result for 
J ~ 850K, corresponding to the smaller gap. 

Troyer et al. calculated x ^i^nd 1/Ti for the ladder by considering the magnon dispersions 
obtained in the limit J2 ^ Ji.i The lowest branch is a single-magnon state which is odd with 

2 



respect to interchange of the two chains {ky = vr). This remains the lowest excitation also 
when J2 = Ji- The smallest gap (A) is at momentum fc^. = vr along the chains. As —>■ 0, 
the one-magnon branch crosses into a multi-magnon continuum. At k^ = the gap is ~ 2 A, 
corresponding to a two-magnon excitation. At low temperatures the thermodynamics of the 
ladder is thus obtained by populating the modes with k^ ^ it, ky = n. The susceptibility 
then has the formi 



X 



y-l/2g-A/T^ (2) 



For T up to ~ A this form is in good agreement with results from exact diagonalizations 
of small systems,0 as well as quantum transfer matrix results.! As mentioned above, the 
agreement with experimental results for SrCu203 is also good, with a A 420K.i 

The NMR spin-lattice relaxation rate is related to the dynamic structure factor S{q, u) 
according toi 

l/Ti = iElAH'5(g,a.-0), (3) 
q 

where Ag- is the nuclear hyperfine form factor. At very low temperatures, the main contri- 
butions to 1/Ti come from momentum transfers ~ 0, = 0, i.e. both the initial and final 
states are on the one-magnon branch at kx ~ tt. Taking into account only these processes, 
Troyer et al. obtained the leading low-temperature formi 

1/Ti ~ I Vope-^/^. (4) 

A behavior close to exponential is seen for SrCu2 03 in the temperature regime lOO/T ^ T < 
At lower temperatures 1/Ti is dominated by impurity effects. The J ~ 1300K 
extracted from fits of (^) to experimental data is markedly different form the J ~ 850K 
obtained from the susceptibility. 

One could certainly argue that SrCu203 is not a perfect ladder system. Most likely, 
Ji is not exactly equal to J2. However, the above theoretical forms only depend on the 
gap, and the disagreement between the gaps from x and 1/Ti cannot be explained by 



Ji 7^ J2 alone. Furthermore, the couphng between the ladders is expected to be weak.i Thus, 
before discarding the single ladder as a good approximation of the system, it is important 
to investigate its behavior in more detail. Whereas the low-temperature form (0) for the 
susceptibility has been verified to be accurate by comparisons with numerical results,!^ the 
form (^) for the spin-lattice relaxation rate has not been tested numerically. At very low T 
it is hard to see why (^) should not apply. However, the temperatures for which the fit to 
the experimental results were made are not very low on the scale set by the gap. It is clear 
that there will be large contributions to 1/Ti from processes with ^ iT,qy = tt between the 
one-magnon branch and the continuum at /c^; ~ if the temperature is high enough for states 
at energies ^ 2A to be populated. These processes are particularly important because the 
ladder has strong short-range antiferromagnetic correlations. The matrix elements entering 
the ^ TT,qy = 7^ processes are therefore much larger than those for ^ 0,qy = 0. Hence, 
although the ~ 0, = contributions are the only ones surviving in the T ^ limit, 
it is quite likely that 1/Ti is actually dominated by other processes at the upper range of 
temperatures considered in the experiments. 

We have calculated 1/Ti using the maximum entropy (ME) methodi to analytically 
continue imaginary time correlation functions obtained by a QMC technique. The spin-echo 
decay rate I/T2G is related to the static susceptibility,^ which can be calculated directly. 
We have used a recently developed QMC method based on stochastic series expansion,0 
which produces results free from systematical errors associated with Trotter based methods. 

The calculations of 1/Ti and I/T2G require knowledge of the Cu nuclear hyperfine inter- 
actions. For the high-Tc cuprates La2_xSr^Cu04 and YBa2Cu306+x the hyperfine couplings 
are well described by the Mila-Rice form,lli with an anisotropic on-site coupling with compo- 
nents y4_|_ and Ay, and an isotropic nearest-neighbor transferred coupling B. Typical values 
reported are B =41kOe//.„ B/A, = 1.2, and B/A, = -0.25.i Knight shift measurements 
on SrCu203 by Ishida et al% indicate that B is much smaller in this compound. Assuming 
that a single ladder picture is appropriate for the spin system as well as for the hyperfine 
couplings, and that the S-couplings have equal strengths along and across the chains, the 



Knight shift results! give the relations 

+ 3Sx = 48kOe//iB, (5a) 
Ay + 35|| = -120kOe//iB, (5b) 

where we have not excluded an anisotropic B. Assuming that the on-site couplings remain 
close to their standard Mila-Rice values, the transferred couplings in SrCu203 are thus Bx_ ~ 
4kOe//iB, and B\\ 15kOe//iB- Considering experimental uncertainties, these estimates are 
probably consistent with 5^ = B\\. In any case, the magnitude of B is much smaller than 
the typical two-dimensional cuprate value. 

Below we present results for x, 1/Ti and I/T2G. For the NMR rates we use the relations 
(§) and several values of the ratios -B_l/A_l and -By /Ay. We consider two choices for the 
exchange J, corresponding to approximately the values obtained before from x and 1/Ti; 
J = 850K and J = 1200K. 

Fig. |1] shows our results for the spin susceptibility along with the experimental results 
by Azuma et a/..i We have used lattices with up to = 2 x 128 spins, which for the 
temperatures considered here is enough for finite-size effects to be negligible. The agreement 
with earlier numerical resultsB'i is good. With a g-factor g = 2 the QMC results for J = 850K 
agree well with the experimental data. For J = 1200K the resulting x fails to reproduce the 
experimental results. Even if g is adjusted it is not possible to obtain a reasonable agreement 
for any sizeable temperature regime. If indeed J is 1200K or larger, this discrepancy is hard 
to explain. 

For extracing l/Ti we have calculated the r-space imaginary time correlation functions 
corresponding to Eq. (^, and continued these numerically to real frequencies using the ME 
technique.i We have obtained the relevant correlation functions to within relative statistical 
errors of 10~^ — 10~^ for systems with up to 2 x 128 spins. Even with this high accuracy 
the continued functions have some uncertainties. At high temperatures the procedures can 
be tested against exact diagonalization results, since the distribution of (5-functions that 
represent the dynamic structure factor of a small system then is dense enough that a small 



broadening produces a smooth function, which can be compared with the results obtained 
with the ME technique. As the temperature is lowered, the number of 5-functions with 
significant weight decreases rapidly. For the largest systems that can be exactly diagonalized 
the presence of many gaps then prohibit meaningful comparisons with ME results, since this 
method cannot resolve structure on that scale. At temperatures where comparisons are 
meaningful, the ME method produces results in good agreement with exact results for a 16 
site Heisenberg chain.El Additional evidence that this is a reliable method for obtaining l/Ti 
stems from work on the two-dimensional Heisenberg model,0 where good agreement with 
experiments on La2Cu04 was found, as well as results for the one-dimensional Heisenberg 
model,lli which exhibit the behavior expected on theoretical grounds. 

For the ladder, results obtained using the ME technique become uncertain at tempera- 
tures where the gap opens up, and the weight for uj ^ relative to the weight for u > A 
decreases rapidly. We beheve that our results are accurate for T > A/2, and become in- 
creasingly inaccurate for lower T. Here we present results for T/J> 0.2. The accuracy of 
the results are probably not higher than tens of percent in the worst cases. Nevertheless, 
they are useful for establishing the general trends. 

Fig. shows results for several values of the ratio Bj_/A^, with relation ( |5aD satisfied. 
Interestingly, for a strictly local interaction {B±/A± = 0) and J = 1200K there is very good 
agreement with the experiment. However, in this case (|5aD gives Aj_ ^ 48kOe//iB; which is 
much higher than one would expect. It is believed that the on-site couplings should be less 
sensitive to details of the structure of a particular material than the transferred couplings, 
and therefore one expects Aj_ ^ 34kOe//iB as in planar cuprates.lll For J = 850K the best 
over-all agreement is obtained with B±/A± ^ 0.1, which gives a reasonable value for A± as 
well. However, the slope of the curve is different from the experimental one. Nevertheless, 
it is interesting to note that the magnitude of 1/Ti agrees with the experimental curve to 
within a factor of 2 in the regime 150K ^ T ^ 300K, with a J = 850K that accounts for 
the susceptibility as well. 

A clear indication that 1/Ti in the regime considered here is not dominated by g ^ 
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processes is that there is a significant decrease in 1/Ti with increasing B±/A±. With 
satisfied, the form factor Aq=Q remains constant, and hence the low-T form (^) predicts a 
1/Ti that does not change with B±/A±. As we argued above, one can expect processes with 
q ^ (tt, 7r) to be important at these temperatures, and the decrease in 1/Ti with increasing 
B^/A±_ is then naturally explained by the decrease in the form factor at g* = (tt, tt). 

Only rough estimates of the behavior of the spin-echo decay rate of the Heisenberg ladder 
have been made.i It is dominated by the indirect nuclear spin-spin interactions induced by 
the coupling to the electronic spin system. Pennington and Slichter derived the formal 
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0.69 



1/2 



E'^KO,^) , (6) 



T2G 

where Jz{xi,X2) is the z-component of the induced interaction between nuclei at Xi and X2: 



Jzixi, X2) = -| E ^(^1 ~ ri)A{x2 - rj)x{i - j), (7) 
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and 0.69 is the natural abundance of Cu isotope. The only non-zero hyperfine couplings 
are ^4(0) = Ay and ^4(1) = B^\. For a system with a gap, the static susceptibility ~ j) = 



dT{Si{T)Sj{0)) decays exponentially with \ fi — fj\ even at T = 0. I/T2G calculated for a 
ladder with 2 x 128 spins at T/ J ^ A is therefore a good approximation to the T = result 
of an infinite system. For this quantity we can thus obtain ground state as well as finite-T 
results. 



Fig. ^ shows results obtained using relation (|5bD and several ratios -By /Ay. If A\\ is 
to remain close to its value in planar cuprates we need -By/^n ~ —0.1, which with J = 
850K indeed gives a quite good agreement with the experimental result. An almost perfect 
agreement is obtained with J = 850K and ^ -0.12. With J = 1200K a slightly 

larger B\\/A\\ is needed to produce an approximate agreement with the experiment, but the 
slope of the numerical curve cannot be reproduced as well as with J = 850K. Note that for 
a strictly local coupling {B\\ = 0) and J = 1200K, which gave a good agreement for 1/Ti 
(Fig. 1^), I/T2G is almost an order of magnitude too small. 

All the above results were obtained with the assumption that the chain coupling Ji is 



equal to the rung coupling J2. It is important to consider also the more general case of non- 
equal couplings. Allowing J2 7^ Ji we find that the best agreement with the susceptibility 
is obtained with J2/ Ji ~ 0.8, and Ji ~ HOOK (this requires a g-factor g 2.1). The 
results for 1/Ti and I/T2G calculated with these parameters are not in significantly better 
agreement with the experiments than those shown in Figs. ^ and 0, however. 

We conclude that the experimentally measured x and I/T2G for SrCu203 can be well 
accounted for by a Heisenberg ladder with J = 850K, and the experimentally determined 
hyperfine couplings. The calculated 1/Ti agrees with the experiment to within a factor of 2. 
The reason for the discrepancies in this quantity could be details of the hyperfine couplings 
not taken into account here, such as possible differences in the transferred couplings B along 
a chain and on a rung. 1/Ti is a direct measure of the low-frequency spin fluctuation spectral 
weight, whereas I/T2G is given by a frequency integral. It is therefore likely that 1/Ti is 
more sensitive than I/T2G to slight deviations from the assumed hyperflne relations in 
the regime where the low frequence spin fluctuation spectral weight drops rapidly. Hence, we 
consider the agreement with the experiment to within a factor 2 reasonable. We propose that 
the reason for the discrepancies reported earlierii for the gaps extracted from x and 1/Ti is 
that contributions to 1/Ti arising from processes with momentum transfer ^ iTjqy = n are 
important at high temperatures. Since only a narrow range of relatively high temperatures 
is accessible experimentally, a flt to the low-T form (^) can give misleading results for A. 

The value of J hence appears to be smaller than the typical values observed in high-T^ 
cuprates. This is puzzling, since the Cu-0 bond structure of the SrCu203 ladders is the same 
as that of the two-dimensional cuprates.0i One possible explanation for the reduced value is 
that J represents an effective coupling once inter-ladder effects are taken into account. The 
weak frustrated ferromagnetic coupling between ladders is expected to enhance the gap,i 
and is therefore not a likely mechanism for reducing the effective J. On the other hand, a 
c-axis coupling reduces the gap, and may be important in SrCu203. However, preliminary 
QMC results for the susceptibility of a stack of weakly coupled ladders with J > lOOOK do 
not compare as favourably with the experiments as the single ladder result with J = 850K 
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shown in Fig |l|. We would also like to point out the possibility that the Madelung potentials 
of SrCu203 and the two-dimensional cuprates might be different. This would alter the size 
of the energy denominators associated with the superexchange interaction. 
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experimental data. This work is supported by the Office of Naval Research under Grant No. 
ONR N00014-93-0495 (A.W.S. and E.D.) and the Department of Energy under Grant No. 
DE-FG03-85ER45197 (D.J.S.). 
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FIGURES 

FIG. 1. QMC results for the spin susceptibility of the Heisenberg ladder with J = 850K and 
J = 1200ir compared with the experimental results for SrCu203. A g-factor g = 2 was used for 
both sets of numerical results. 

FIG. 2. The spin-lattice relaxation rate calculated using QMC and ME compared with the 
experimental results by Ishida et a/.i (thick solid curves). The upper and lower panels show results 
for J = 850K and J = 1200K, respectively. The hyperfine couplings used satisfy (|5a|). The ratios 
B±/A± are (open circles) 0.05 (solid circles), 0.10 (open squares), and 0.20 (solid squares). 

FIG. 3. QMC results for the spin-echo decay rate compared with the experimental results by 
Ishida et a/.i (solid curves). The upper and lower panels show results for J = 850K and J = 1200K, 
respectively. The hyperfine couplings used satisfy the relation (|5b|). The ratios -By /Ay are (open 
circles) —0.05 (solid circles), —0.10 (open squares), and —0.15 (solid squares). The dashed curve 
is the best fit for J = 850K, with B»/A» = -0.12. 
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Fig. 3, A. W. Sandvik et al. 



